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Abstract 

Through an explicit calculation of massive scalar particle (e.g., Higgs boson) production in high 
energy pA collisions up to one- loop order, we demonstrate that the Sudakov- type logarithms in 
hard processes in small-x saturation formalism can be systematically separated from the typical 
small-x logarithms. The generic feature of the Sudakov logarithms and all order resummation is 
derived. We further comment on the phenomenological implications and extension to other hard 
processes in the small-x calculations. 
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1. Introduction. An important application of the perturbative quantum chromo dynamics 
(QCD) is the resummation. In high energy hadronic processes involving large separate scales, 
resummation is not only necessary to make reliable predictions, but also crucial to extract 
the fundamental properties of the strong interaction theory, such as the strong coupling 
constant. In practice, the resummation formalism has been applied to a wide range of 
physics processes. One of the examples is the resummation of the Sudakov-type double 
logarithms [H, H]- The double logarithms appear in, for example, the transverse momentum 
spectrum of a hard process, where each order of perturbative correction is accompanied by 
a large double logarithmic term of \ii'{Q'^ /k\) with Q the large momentum scale and k± the 
transverse momentum. In low transverse momentum region where most of the production 
events sit, the QCD resummation has to be performed. This resummation is often referred 
as the transverse momentum resummation. 

Meanwhile, there is also small-x resummation which is equally important, in particular, in 
the area of the large hadron collider (LHC). Small-x resummation is governed by the BFKL 
evolution jsj. Because of high gluon density in nucleon/nucleus, the non-linear term in the 
evolution plays a very important role at small- a; 0-[6|, which leads to the BK-JIMWLK 
evolution |7H9|. As a result, the gluon saturation with a characteristic scale becomes 
inevitable at small-x. Seeking for the signal of the gluon saturation phenomena and studying 
the associated dynamics has been one of the most important motives for the high energy 
nucleon-nucleus experiments at RHIC and the LHC, and for the planed electron-ion colliders. 
Great efforts have been made from both experiment and theory sides [lo|. Among them, 
the hard processes involving a large momentum scale have been emphasized recently as 
important probes for the saturation phenomena This is because these hard processes 
can directly measure the transverse momentum dependence of the gluon distributions, whose 
behavior manifest the saturation phenomena js], |6[. However, these hard processes also 
impose a question: are there Sudakov double logarithms as well? If yes, how to resum these 
double logarithms consistently in the small-x saturation formalism? To our knowledge, there 
has been no theoretical investigation on this topic. 

In this paper, we will study, for the first time, the Sudakov-type double logarithms in 
the small-x formalism. As an important first step, we take the massive (with mass Q) color 
neutral scalar particle production in pA collisions as an example. The scalar particle is 
directly coupled to the gluon through an effective lagrangian £e// = ~\g<j>^F^u^'^^'' ^ where 
$ represents the scalar field, F"-^^ is the field tensor for the gluon field with the associated 
color index a. This effective lagrangian has been used to calculate the Standard Model 



Higgs boson production in hadronic collisions [12[, and to study the gluon saturation in 
nucleus [13] as well. We extend the previous calculations to the scalar particle production in 
pA collisions in the small-x factorization formalism up to one-loop order. At this order, we 
will be able to identify the Sudakov double logarithms, which is absent in the leading order 
evaluation. As a result, all order resummation can be performed consistently with the small- 
X evolution. One of the important steps to achieve this is to separate the soft gluon radiation 
(which contributes to the Sudakov logs) from those contributing to the small-x evolution. 
Our explicit one-loop calculations will demonstrate that we can consistently resum Sudakov 
double logarithms and the BFKL (BK-JIMWLK) evolution at the same time. This will 
provide important guidelines for further developments in other hard processes which are 
crucial to study the saturation phenomena in the small-x physics. 

The rest of the paper is organized as follows. In Sec. II, We will present the calculation 
of scalar particle production in pA collisions in the small-x formalism up to one-loop order. 
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FIG. 1. Examples of one- loop diagrams from the real gluon radiation (L) and virtual gluon radiation 
(R) for the scalar particle production in pA collisions in the saturation formalism with the multiple 
interactions with the nucleus target taken into account. Here the vertical gluon represents arbitrary 
number of the WW small-x gluons which are summed into Wilson lines. 

The Sudakov double logarithms are identified. All order resummation is derived in Sec. III. 
We summarize our results and present further discussions in Sec. IV. 

2. Massive Scalar Production in pA Collisions at One-loop Order. We follow the high 
energy small-x factorization formalism to calculate the massive scalar particle production 
in pA collisions which implies that the center of mass energy S ^ while ^ k"^. 
Similar to those calculated in Ref. p^], the gluon distribution from the nucleon is taken to 
be coUinear, therefore the transverse-momentum distribution of the produced scalar particle 
reflects the transverse-momentum dependence of the gluon distribution in the nucleus. Due 
to the high gluon density inside the target nucleus, the multiple interactions must be taken 
into account. 

The lea ding order contribution can be easily formulated following the small- a; factorization 



formalism 13 



c/a(LO) f d^x^d^x' 



.„y ^^e--.j,wsr-fe.x'j , (1) 

where r± = x± — x'j_, ao = g"^/ g^32{l — e) with e = — (D — 4)/2 the dimensional regulation 
parameter, y and k± are the rapidity and transverse momentum for the scalar particle, gp{x) 
denotes the integrated gluon distribution from the nucleon, and represents the uninte- 
grated gluon distribution function from the nucleus. Here, x ~ Qe^ / ^/S and Xa ~ Qe^^/ 
are the momentum fractions of incoming gluons from the proton projectile and target nu- 
cleus, respectively. Because of the colorless nature of the scalar particle, there are only the 
initial state interactions for the multiple gluon exchange with the nucleus target. Therefore, 
the associated un-integrated gluon distribution is the so-called Weizsacker- Williams (WW) 



gluon distribution in the small- a; formalism [11 



S^'^{x^,y^) = - (Tr [diU{x^)U\y^)d^^U{y^)U\x^)])^ , (2) 

where Y represents the rapidity of the gluon from the nucleus Y ~ ln(l/xa), and sum over 
the transverse index /3 is implicit. The so-called linearly polarized gluon distribution can also 
be taken into account where different projection of the transverse index is performed (l5| . 

The Wilson line U is defined as U{xj_) = V exp (^—ig dx~A~^{x~ ,x±)^. 

At one-loop order, real and virtual diagrams yield higher order corrections which contain 
various divergences. In addition to the rapidity divergence associated with the WW gluon 
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distribution from the nucleus, there is also the collinear divergence for the integrated gluon 
distribution from the nucleon. In our calculations, the dimensional regulation (D = 4 — 2e) 
is used to formulate the collinear divergence, whereas the direct subtraction is applied for 
the rapidity divergence. In Fig. [H we plot the typical diagrams from the real and virtual 
gluon radiations at one-loop order . The scattering amplitude from the real diagrams can 
be written as, 



+ 



av 



k'l + ej kl + e})^ ^ 2 \kl + e) k'f + 



TB{kg±) , (3) 



where k± and k2i_ represent the transverse momenta for the final state scalar particle and 
the radiated gluon, respectively, kgx^ = kgi±_ + kg2± = k± + k2±, = {1 — z)Q'^ with z the 
momentum fraction of the incoming gluon carried by the scalar particle. For convenience, 
we further define ki± = k± — kgi±, k'j_ = k± — zkg±, and z/ and a are the polarization 
vector indices for the incoming and outgoing gluons, for which we have chosen the physical 
polarizations. and F^ are defined as 

rl(fc3ix,fc,2±) = / rfVc^'xae^'^^^-^^+^'^^^^-'^Tr [f/t(x2)T^f/(x2)[^9^^f/^(xi)f/(xi),r^l ,(4) 



TBik 
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where a and b are color indices for the incoming and outgoing gluons, respectively. Clearly, 
the amplitude squared from the above expressions will depend on multi-gluon correlation 
functions (beyond the WW-gluon distribution) from the nucleus, as this is the common 
feature in the high order calculations in the small-x formalism [l^ . However, in the k± <^ Q 
limit, these correlation functions is either reduced to the WW-gluon distribution, or absorbed 
into the evolution of the WW-gluon distribution. To evaluate the contribution from the real 
gluon radiation, we integrate out the phase space of the radiated gluon (/c2)- To simplify the 
calculation, we perform the power expansion of the amplitude squared in terms of k±/Q, 
and only keep the leading power contributions. In the power counting analysis, we find that 
the phase space integral contains three important contributions: (1) soft gluon radiation 
A;^ ~ A;^ ~ which eventually leads to the Sudakov logarithms; (2) collinear gluon 
contribution as respect to the incoming nucleon projectile; (3) collinear gluon contribution 
as respect to the target nucleus. The soft gluon contribution can be easily obtained in the 
limit of ^ k'j_, which results into ln(Q^/fc|_|_). When Fourier transformed 

into the impact parameter space, this term leads to a soft divergence in terms of in 
dimension regulation. The soft divergence will be cancelled by the relevant virtual diagrams. 
The last contribution contains the rapidity divergence and gives rise to the evolution of the 
WW gluon distribution. 

The evaluation of the virtual diagrams leads to the following contributions. 



where e'j = —z{l — z)Q'^, q'± = (1± + kgi± — zk± with k± = kgi±_ + fcc,2±, /3o = 11/12 — Nf /18 
with Nf number of flavors, euv and em represent the uhra-violet and infrared divergences 
in the loop diagrams, respectively, and Tc and Yd are defined as 

Tc{kgi±,kg2±) = J d^Xid^X2e''o'^-'''+''o^^-''' 

X {Tr [f/(xi)T"f/t(x2)] Tr [U^ {x^)Uix2)] - (xi ^ 0:2)} , (7) 
Tfj{k^) = j d\e'^o^-''^Ti [T,i9^[/t(xx)f/(xx)] • (8) 

In Eq. (E]), the ultra-violet (UV) divergences cancel out between the two contributions in the 
first term, whereas the remaining UV- divergence in the second term is normally interpreted 
as the charge renormalization, which appears in the form of {as/7i)Ncf3o (— l/e^/y + ln(Q^//i^)). 
Furthermore, there is also IR divergence, which can be absorbed into the renormalization for 
the incoming gluon distribution from the nucleon. In addition, there are rapidity divergence 
and soft divergence in the first term of Eq. ([6]) . The soft divergence is regulated by dimension 
regularization, and is found to proportional to 1/e^. This divergence, as mentioned above, 
will cancel out that from the real diagrams. 

After canceling out the soft divergences from the real and virtual diagrams, we are left 
with the rapidity divergence and coUinear divergence, 

xgp{x) j yKoMMX ® + ^^^S^^ {x±,y±)Vg/g <^ xgp{x) , (9) 

where Kdmmx is the BK-type of evolution kernel for the WW gluon distribution, and Vg/g = 
"t" i + ~ ^) + /5o'^(l ~ z) the collinear evolution kernel for the incoming gluon 
from the nucleon. Our result for Kdmmx is consistent with a recent calculation for the 
evolution directly from the JIMWLK evolution [16]. After subtracting these divergences 
by the renormalizations of the associated distributions, we obtain the final result for the 
differential cross section as 

^^(LO+NLO) r ^2^ ,2 



dyd'^kj 



f d X±d ifcx-r^ CWW ( ^' 1,2 _ 



1 + 

IT 



Cq 2 \ Cq 



(10) 



where cq = 2e~'^'^ with 'Je the Euler constant. We have also introduced the scale for the 
integrated gluon distribution /i^ = Cg/r^ and the rapidity for the WW gluon y = In 1/xa to 
simplify the above expression. 

Comparin g to the differential cross section calculation for inclusive hadron production in 
pA collisions [lj|, we find that they share the same structure. In particular, the integrated 
parton distribution is more convenient to be set at the scale of co/r±. The major difference 
is that, the above formula is valid in the limit of Q S> fc^ where we have to resum the 
Sudakov logarithms. That is also the reason that the un-integrated gluon distribution from 
the nucleus side does not depend on more complicated structure of the Wilson lines. If we 
calculate the differential cross section in the kinematic region of fc^ ~ Q, the multi-gluon 
correlation functions (beyond the WW-gluon distribution) from the nucleus will contribute. 



similar to what was found in Ref. [IJ]. Those contributions, however, is power suppressed 
in the limit oi k± Q. 
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We have also performed this calculation in coordinate space with cutoffs putting into 
integrals instead of the use of the dimensional regularization, and obtained the same results. 
This implies that there is no scheme dependence in the Sudakov logarthms, at least to the 
level we are concerned. 

3. All Order Resummation. Our calculations at one-loop order in the above demonstrate 
that the soft gluon radiation is well separated from the collinear gluon radiations. In par- 
ticular, the soft gluon radiation comes from the initial state radiation, and in the Q ^ k± 
limit, the relevant gluon distribution from the nucleus takes the form of the WW gluon 
distribution. Higher order soft gluon radiations will follow the same form. This can be 
understood that the Sudakov double logarithms come from soft gluon radiation associated 
with the hard probe. To resum these large logarithms, we follow the Collins-Soper-Sterman 
procedure j2|. In particular, we can write down an evolution equation respect to the hard 
scale Q^. By solving the differential equation, we can resum the differential cross section in 
terms. 



^^(resum) 



(27r)2 

xxgp{x,iJ,'^ = cl/rl) 



1 + — — iVc 
TT 2 



where the Sudakov form factor contains all order resummation 



s 



sud 



(11) 
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where Ci^2 are parameters in order of 1. The hard coefficients A and B can be calculated 

oo 

perturbatively: A = ^ A*^*^ iirT ■ From the explicit results for the one-loop calculations, 

1=1 ^ 

we find that they are 



(13) 



where we have chosen the so-called canonical variables for Ci = Cq and C2 = 1. The 
above formulas are the final results for the soft gluon resummation in massive scalar particle 
production in pA collisions. The result is valid in the limit of k± -C Q, and we have applied 
the small- a; factorization where higher order in l/ln(l/a;a) has been neglected as well. As 
compared to Ref. which takes the as correction of g^i, into account, we find the above 
coefficients are consistent with the known resultsj3] except for B^^^ which misses a factor of 
2. This difference is due to the convention in the saturation formalism in which we do not 
include the virtual gluon and quark loops contribution for the vertical WW small-x gluons 
as shown in Fig. [H 

We would like to emphasize a number of important features in the above derivation. First, 
the collinear divergence and the soft divergence is well separated. This follows the transverse 
momentum resummation in the collinear factorization. Second and most importantly, the 
rapidity divergence from the the un-integrated gluon distribution of the nucleus is also well 
separated from the soft gluon radiation. This is because the rapidity divergence comes 
form the collinear gluon radiation parallel to the nuclei momentum, whereas the Sudakov 
logarithms come from soft gluon region. 

Third, as the above features are common features in the small-x saturation formalism, 
we expect our resummation results can be extended to other hard processes. An immediate 
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application is the color neutral particle production in pA collisions, such as the two-photon 
production [l^ and heavy quarkonium production [20]. There are only initial state interac- 
tions in these processes, and the Sudakov resummation will be the same as what we have 
derived in this paper for the massive scalar particle production. We expect the similar resum- 
mation formula. Another important extension is the dijet correlation in pA collisions [TT], 
which has attracted great attentions in recent years. Of course, because the final state of 
dijet production carries color, we have to take into account the final state interaction effects 
in the Sudakov resummation. Therefore, the resummation formula will be different from the 
scalar particle production studied in this paper. However, the generic form will remain the 
same. We plan to address this calculation in a future publication. 

4. Discussions and Conclusion. In summary, we have demonstrated the Sudakov double 
logarithms in the small-x calculations of hard processes in pA collisions. All order resumma- 
tion formalism has been derived. This result will have potential application for other hard 
processes in pA collisions, from which we hope to investigate the saturation physics. 

In addition, the technique developed in this paper shall also be relevant for the hard 
processes in hot/dense medium, such as the quark-gluon plasma. In particular, in the hard 
processes of jet penetrating through the medium, we expect the similar formalism for the 
Sudakov resummation effects. The direct consequence is that the transverse momentum 
broadening in the hard processes (with hard momentum scale much larger than the mo- 
mentum scale in the medium) will be dominated by the Sudakov logarithms which is the 
same as that in the vacuum. This may lead to a natural explanation for the azimuthal angle 
correlation of dijet production in heavy ion collisions which was found the same as that in 
pp collisions in the large angle region j2l|, where Sudakov effects dominate. 
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DE-AC02-05CH11231 and DOE OJI grant No. DE - SC0002145. We thank J. W. Qiu for 
comments and discussions. 
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